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Periodic in time spin-Hall current and spin polarization induced by a dc electric bias has been cal- 
culated in a superconductor-normal 2DEG-superconductor (SNS) Josephson junction. We assumed 
that the band energies of electrons in the normal system are splitted due to Rashba spin-orbit cou- 
pling. The transport parameters have been calculated within the diffusion approximation and using 
perturbation expansion over a small SN contact transparency. We found out that in contrast to the 
stationary Josephson effect, the spin-Hall current does not turn to zero. Besides a direct proximity 
effect caused by Cooper pair's transition into a triplet state, the spin current and polarization are 
also driven by a periodic electric field associated with the charge imbalance. 

PACS numbers: 72.25.Dc, 71.70.Ej, 73.40.Lq 
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I. INTRODUCTION 

The spin-Hall effect (SHE) is a fundamental physi- 
cal phenomenon where the spin-orbit interaction (SOI) 
shows up in electron transport on a macroscopic level. 
The interplay of spin precession caused by SOI and elec- 
tron acceleration in the electric field gives rise to a flux of 
the out-of-plane spin polarization flowing perpendicular 
to the electric current. Although this effect was predicted 
long time agC''^ , it has been observed experimentally only 
recently in semiconductors^''^ and metals^. The nature 
of this phenomenon is now well understood and stud- 
ied for various systems (for a review see (H) . The spin- 
Hall effect is being considered as a tool for manipulating 
electron spins in perspective spintronic applications. On 
the other hand, the spin polarization accumulated due to 
SHE is subject to dissipative processes of spin relaxation 
and diffusion. From this point of view, it is interesting 
to consider SHE in superconductors, as well as in SNS 
junctions, where the N-region is represented by a normal 
electron system with a strong enough spin-orbit coupling. 
Also, in such systems SHE provides an opportunity for 
a direct coupling of spin degrees of freedom to supercon- 
ducting quibits^''' 

An important distinction of spin-Hall effects in super- 
conducting and normal systems is that in the latter case 
this effect is determined by spin dynamics of single par- 
ticles, while in the former case major role is played by 
interference of triplet and singlet Cooper pairsS. The 
triplet correlations, in their turn, are induced in the con- 
densate wave function by SOI^, and shows up as admix- 
ture to the singlet state. Therefore, the spin current and 
spin accumulation caused by SHE are determined by a 
coherent macroscopic state and do not dissipate. SHE 
has been considered in bulk superconductors and SNS 
tunneling contacts^. Also in such contacts an effect recip- 
rocal to SHE was recently studied^ ^. These studies have 
been restricted to the stationary transport. In the case 
of SNS junctions this means that the Josephson electric 



current is driven by a phase difference of superconduct- 
ing order parameters of two superconducting electrodes. 
If these electrodes have different electric potentials, this 
current will periodically vary in time. One would expect 
that, due to such a time dependence, the non-zero spin- 
Hall current will be induced, while it is forbidden in the 
stationary case by the time inversion symmetryi^ One 
more nonstationary effect is associated with an electric 
fleld caused by a dynamic electron-hole charge imbalance 
within the normal layer. Such a periodic field will drive 
a flux of the spin polarization, in a way quite similar to 
conventional SHE in normal systems. 

The goal of the present study is to extend the theory of 
RefH to the case of the non-stationary Josephson effect. 
We will calculate the spin-Hall current and spin polariza- 
tion created by a combined effects of the Josephson tun- 
neling and SOI. It will be assumed that two singlet super- 
conducting electrodes are under the dc electric voltage V. 
The normal layer is contacted to them through the low 




FIG. f : (Colour online) An SNS Josephson junction. The 
ac Josephson electric current flows in a;-direction between 
two superconducting electrodes (blue) through normal 2DEG 
(gray). Black layers show tunneling barriers. In the pres- 
ence of Rashba spin-orbit coupling the ac spin current of 
z-polarized electrons in y-direction is induced, with zero z- 
oriented spin density. A finite spin polarization in y-direction 
(not shown) is also induced . 
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transparency tunneling barriers, as shown in Fig.l. This 
layer is taken thin enough, so that electrons are restricted 
to 2D motion, as, for example in a semiconductor quan- 
tum well. The spin-orbit coupling in the N-layer is given 
by the Rashba interaction^^ and wc will ignore the spin- 
orbit effects in scattering of electrons from impurities. 
At the same time the spin-independent scattering will be 
taken into account within the Born approximation. The 
particle's mean free path I will be assumed smaller than 
all relevant parameters of length dimension, except the 
Fermi wavelength kp, which in the semiclassical approx- 
imation is much smaller than I. Therefore, the electron 
transport within the N-layer is dominated by diffusion. 

The outline of this paper is as follows. The general 
expressions for the spin-Hall current and spin density are 
derived in Sec. II. In Sec. Ill some numerical results are 
presented and discussed. Finally, Appendix A presents 
some details of analytical calculations within the Keldysh 
formalizm. 



II. BASIC EQUATIONS 

Since we will focus on the basic characteristics of SHE, 
the simplest approach will be employed within the lowest- 
order perturbation theory with respect to transmission 
coefficients of interface barriers. It should be noted, how- 
ever, that higher-order corrections to the electric Joseph- 
son current are not always small^^, in particular in the 
range of temperatures T larger than the Thouless energy 
Erh — D/L^, where D is the diffusion constant and L is 
the distance between contacts. We will assume the tem- 
perature or/and the transmission coefficient low enough 
to avoid such a situation. The tunneling is presented by 
the perturbation Hamiltonian 

Hint ^ ^ tk,k'a\,^TzCkcr +h.C., (1) 
k.k\a 

where aka and Cka are electron destruction operators in 
the superconductor and N-layer, respectively, with k de- 
noting the wavevector and a the spin projection of the 
particles. The transmission coefficient will be assumed a 
slow varying function of wavevectors in the vicinity of the 
Fermi surface. Hamiltonian ([T|) is written in the Nambu 
representation, where destruction operators are defined 
as 



and Ti, T2, T3 are the Pauli matrices in the Nambu space. 
In its turn, the unperturbed Hamiltonian of the normal 
layer has the form 

k,a,<j' \ / 

Uk,k'c\,„T'iCka + He, (3) 



where cr = [a^.a^.a^] is the Pauli spin vector. The 
Rashba spin-orbit field hk, which is the odd function of 
k, is given by^^ fv^ = aky, fv^ = —ak^. The random im- 
purity scattering potential is represented by its matrix 
elements Uk,k'- This scattering determines the elastic 
mean free time r of electrons. For a short-range scatter- 
ing it is given by l/r = 2TrNF{\Uk.k' P)imp = 27rn,7VF|C/p 
, where Np is the state density at the Fermi level, Ui is 
the impurity concentration and the subscript "imp" de- 
notes averaging over impurity positions. The Hamilto- 
nian He represents the Coulomb interaction of electrons. 
It will be treated within the random phase approxima- 
tion to take into account screening effects associated with 
the dynamic charge imbalance, while its contribution to 
the self-energy and electron-electron correlations will be 
ignored. 

We assume that SNS contact is unbounded in the y- 
direction. Hence, the Josephson current is in the x- 
direction, as shown in Fig. 1, and the spin-Hall current 
polarized parallel to the z-axis flows in the y-direction 
and depends on the x-coordinate. In the framework of 
the Keldysh formalism^ it can be written as 

Ux,t) = ^(V^, - V^)Tr[a-(Gfi(i,r;t,r'))imp]|r^r', 

(4) 

where G'^(t, r; t, r') is a nondiagonal (Keldysh) compo- 
nent of the Green function, which is a 2 x 2 matrix in the 
spin space, with subscript 11 denoting the corresponding 
projection in the Nambu space. Besides the spin cur- 
rent we will calculate also the spin polarization along 
the y-axis. In normal systems such a spin polarization 
is due to the electric spin orientation^^. It is usually 
associated with SHE and takes place also at stationary 
Josephson tunneling conditions*. This polarization can 
be expressed as 

Syix,t) = -^Tr[a^(Gfi(t,r;t,r))i^p] (5) 

Expressions Q and (O have to be expanded up to the 
4-th order with respect to the transmittance. The rele- 
vant diagrams are shown in Fig. 2. Comparing them to a 
diagram representation of the charge Josephson current^^ 
one can see that in the latter case diagrams are much 
simpler. That is because conservation of the Josephson 
current allows to reduce its calculation to calculation of 
the charge time derivative in one of the superconducting 
terminals. Unlike the Josephson current, the spin cur- 
rent does not conserve and such a simplification is not 
possible. An additional problem is caused by time de- 
pendence of the charge transport. It results in an electric 
potential inside the N-layer. Therefore, one should take 
into account diagrams which explicitly take into account 
the Coulomb screening there. These diagrams are shown 
in Fig. 3. The calculations below will be restricted to 
the low temperature T and small enough DC voltage V, 
both much less than the superconducting gap A. In this 
regime the electron transport through the contact will 
be dominated by tunneling of Cooper pairs between two 
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FIG. 2: (Color online) Examples of Feynman diagrams for 
calculation of the spin current and spin polarization. The op- 
erator O denotes the spin current, or spin density operators. 
Grey boxes are diffusion propagators. Boxes denoted by "K" 
and "D" relate to Cooperon and diffuson, respectively, "r" 
and "a" stand for retarded and advanced Green functions. 
Slashes denote the tunneling amplitude 



superconducting electrodes, while a dissipative transport 
due to electronic excitations will be exponentially sup- 
pressed. Henceforth, the Green functions of supercon- 
ductor terminals are represented in Fig. 2 and Fig. 3 
by corresponding anomalous functions. More details on 
calculation of Feynman diagrams in Figs. 2 and 3 can be 
found in Appendix A. 

The main building blocks of diagrams in Fig. 2 and 
Fig. 3 are unperturbed equilibrium Green functions aver- 
aged over impurity positions. For the normal layer these 
functions are determined by Hamiltonian ([3]) and are rep- 
resented by their retarded (r) , advanced (a) and Keldysh 
components 



where Eu = {k"^ 12771*) - ^ and F = 1/2t, 



(6) 



G°^(a;,k) = tanh 



2fcsr 



G°'-(w,k)-GO''(w,k) . (7) 



Important entries in Fig. 2 are the propagators D and K 
given by 



Daf3i,f_t{u!l — UJ2) = 

n.|;7|2(G;^^ii(r,r',c.i)G^^;n(r', 



,W2))i 



(8) 



FIG. 3: (Color online) Examples of Feynman diagrams for 
calculation of Coulomb screening effects on the spin current 
and spin polarization (see Fig. 2 for details). Dashed line 
denotes the screened Coulomb potential. 
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The conjugated functions K^' are defined by Eq. Q with 
interchanged Nambu subscripts 11 and 22. Within the 
semiclassical approximation D and K can be represented 
by the ladder series.— At small frequencies and large 
|r — r'l :S> I they satisfy a diffusion equation and are 
called " diffuson" and " Cooperon" , respectively Due to 
the time inversion symmetry these correlators are not 
independent. They can be expressed via each other. 
Depending on a combination of spin components, the 
diffuson and Cooperon describe either spin, or particle 
(charge) diffusion. Therefore, it is convenient to expand 
them in terms of Pauli matrices, according to 



E>al3v^ 



1 
2 



(10) 



where i,j = x,y,z,0 and a" denotes the 2x2 unity ma- 
trix. Here and below a summation is assumed over the 
vector, or spinor indexes entering twice into an expres- 
sion. Cooperon components can be represented in a way 
similar to (|10l) . The tensor components Dij have a clear 
physical meaning. Namely, Dqq relates to the particle 
diffusion, while various components with i, j = x,y,z are 
associated with the spin diffusion. Mixed terms, for ex- 
ample DiQ, are generally not zero in the presence of SOI. 
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As follows from definitions ([2]) and ([9]), z-components of 
K are related to diffusion of singlet Cooper pairs, be- 
cause they involve antisymmetric combinations of "up" 
and "down" spins in the particle-particle scattering pro- 
cess associated with the propagator K. Other compo- 
nents, 0,x,y are related to triplet Cooperons. It is easy 
to see that the the 0-term gives a triplet with a spin pro- 
jection onto the z-axis, while x and y components are var- 
ious combinations of ±1 triplets. A singlet-triplet mix- 
ing is associated with nondiagonal correlators K^^n, K^z, 
where m = 0, ±1. As it will become clear below, the mix- 
ing terms are proportional to the small parameter ajvp^ 
where vp \^ the Fermi velocity. Only linear in this pa- 
rameter terms will be taken into account in the following 
calculations of the spin-Hall current and spin polariza- 
tion. 

Since K always enters together with the Green func- 
tions of superconducting terminals, it is convenient to 
introduce the pairing function 

(11) 

where are determined by the anomalous supercon- 
ductor Green functions G'']^2°^ ^e\\ as by geometry of 
contacts and their transmittance. Similarly, the conju- 
gated functions i^t''(a) g^j.g defined through and . 
We assume that the SNS junction is symmetric, with the 
electric potentials ±y/2 applied to the left and right elec- 
trodes, respectively. These potentials result in the time 
dependent factors exp(±if7(i + t')/2) in superconductor 

functions G'i2(i,i') and G§i(i,t')> where VI = ±eV . In 
this case, / can be written in the form 



/ 



A [c{x - XL)Sn,eV + c{x - XR)Sn-ev] 



(12) 



where a can be expressed^ through the resistance Rb 
of the SN interface, as a = l/Ae^NpRb and c{x) is de- 
termined by a profile of the contact. For simplicity, as- 
suming that the distance between contacts L — xr ~ xl 
is much larger than the contact length, c{x) will be ap- 
proximated by the delta-function. Since for our choice 
of parameters A ^ w, fi, the retarded and advanced 
anomalous functions / coincide, we will skip the labels 
r, a below. The functions are obtained from Eq. 
by the substitution L R. 



Let us introduce the vertex function 



k 

G°;ii(w + f^, k + ci)Dp^,^in, q)al, . (13) 

Further, according to the diagram representation in Fig. 
2 and Fig. 3 the spin-Hall current can be written as 



Js{x,t) 



E 



q,q' ,n=±eV 



did 

2^ 



[Jl + J2- M2n, Q)Vq'-^^ {Jlscr + J2scr)] , (14) 



where Q = q -\- q' , while Ji and J2 are given by dia- 
grams (a) and (b) in Fig. 2, respectively. Other terms 
in the integrand take into account Coulomb screening, 
as depicted in Fig. 3. Vq denotes the screened Coulomb 
potential 



Va 



e{2n,q) 



(15) 



At small 51 ^ Ep and q ^ kp the dielectric function 
e(2f2, q) is represented by the hydrodynamic expression 
(see e.g. Ref. [H). 

e{2n,q) = l + Vq^Dq^Doo{2n,q), (16) 



where the diffusion propagator DQo{uj,q) is given by 
Eq. p8|. Taking into account that the two-dimensional 
Coulomb interaction Vq — 27r/eo(7, it follows from Eqs. 
(fTB)l and P5|) that at small il and q the second term in 
Eq. p6l) dominates. Retaining in e only this term we 
arrive to 



VqT 



" T + VqNFDq^Doai2n,q) NpDq^Doo{2n,q) 

(17) 

This approximation corresponds to a complete screening 
of charge within the length scale much larger than the 
screening length. 

Using the above definitions of Green functions and cor- 
relators, various terms in the integrand of Eq. can 
be expressed in the form 



where Fi^q are spacial Fourier transforms of ^-coordinate 



dependent functions Fi defined as 
F, = l=Tr[a'F] 



(19) 
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and the coefRcients buj are given by 



(20) 



Each of the symbols a,b,c take the values r, or a, while 
the subscripts ran through 0,x,y,z. Introducing 

also the coefficients 



with A defined by the expression 

A?i^-^(q, q') = E G'^Ulu -n,k- q') 



abc 



(22) 



aiG'l'liuj + n,k + q)(7,G;j^(a;, k)^ 



(21) other terms in Eq.(fT4|) are expressed as 
I 



J2 = - ji (217, q + q')c?J'J'F[„ (uj + -)f\(Lu--)( tanh — tanh 



\M2n, q + q')cr™/,(^ + ^)^!^, i. - §) (tanh ^ - tanh ^ 



(23) 



^i.c. = Ic^FTiu. + ^)F^:, {u. - ?) tanh ^ - ^c-^;-^ + ^)f]% - ^) tanh " + " 



2kBT 



2kBT 



(24) 



J2scr = ^I)oo(217, q + q')C:F/;,(a. + ^)fl, - ^) (tanh ^ - tanh 

ii5oo(2f7, q + q')cS™/,(^ + ^)^!:' - (tanh - tanh ^) , (25) 



The spin density given by Eq. ([5]) is calculated in a way 
similar to the spin current. The same equation as Eq. 

can be used with the following changes: in Eg. p^ 
the factors Cyij from Eq. ([22|) should be used instead of 
bzij] in Eq. ([23|) one should substitute 2Dyi{2^l,q) in- 
stead of j;(2ri,g). Also, in Eq. ([T4|) the vertex jo{2Q,q) 
should be substituted for 2Dyo{2Q,q). 

Before proceeding with further calculation of the spin- 
Hall current and spin density, it is useful to discuss the 
physical meaning of Eqs. ([TH]), (gS]) and Ji 
gives a "bulk" contribution to the spin current (spin 
density). It is determined by diffusion of Cooper pairs 
from the left and right superconducting leads. Since for 
a chosen range of parameters the distances from the leads 
are much larger than the coherence length in the normal 
metal ^D/A, where D is the diffusion constant, the pen- 
etration depth of Cooper pairs into the normal metal is 
determined by the diffusion length during the time much 
larger than A~^. Thus, the characteristic diffusion time 
of singlet pairs is of the order of min[(fcBT)~^, {eV)~^], 
while in the case of triplets the spin relaxation time comes 
into play, if it is shorter than the diffusion time of singlets. 
Therefore, if the spin relaxation time is shorter than the 
Thouless time E^^ — L^/D, the triplet components Fi 
(i = 0, X, y) of the pairing function will be localized rel- 
atively close to the leads. At the same time, at Eth ^ 



max[fcsT, eV] the singlets F^ can propagate through en- 
tire junction. As it will be shown in the next section, Ji 
is represented by a combination of a pure singlet term of 
the form F^F^ and singlet-triplet interference contribu- 
tions, like FzE^. The former can penetrate over large 
distances, independent on the magnitude of the spin- 
relaxation rate associated with the spin-orbit coupling. 
Therefore, at low enough T and V the spin current rep- 
resented by Ji in Eq. (|14p can be observed far from the 
contacts. At the same time, the spatial distribution of 
the current given by J2 is determined by the spin den- 
sity created by one-particle spin diffusion near one of the 
contacts. This diffusion is represented by the diffusion 
propagator D in Fig. 2b. The diffusion in this case is 
restricted by the spin relaxation length. If this length is 
less than L the corresponding spin current (spin density) 
will be distributed relatively close to contacts. So, it is 
of the "surface" type. The remaining screening terms 
in Eq. (|14p are determined by the long-range Coulomb 
interaction. Therefore, their contribution will be of the 
"bulk" type. 
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III. RESULTS AND DISCUSSION 

In this section we will calculate the functions and co- 
efficients entering into general expressions Eqs. (ITil) - (P5)) 
and present numerical results for spin transport param- 
eters. 

The matrix Dij can be found from the diffusion equa- 
tion. In the SHE regime this equation has been derived in 
a number of works The Cooperon Kij can, in its turn, 
be expressed through Dij. Since the latter depends only 
on the a:-coordinate, the diffusion drift of particles occurs 
in the x-direction. Hence, the effective "magnetic" field 
hk induced by the Rashba interaction is directed paral- 
lel to the y-axis. Therefore, electron spins precess in the 
zx plane. This means that one has coupled equations for 
Dzj and D^j , while Dj,^ -components are decoupled from 
them. They, however, stay coupled to the "charge" dif- 
fuson Doj through the weak spin-charge coupling . For 
example, the mixed function DyQ{uj, x) satisfies the equa- 
tion 

- lojDyo - D\/lDyo + TsoDyo - 2rxV,i^oo - , (26) 

where x — ^^oTso/'iF'^ and the D'ykonov-Perel'^spin re- 
laxation rate is Tso — a^kp/T. The spin-charge coupling 
X will be taken into account in the lowest order pertur- 
bation expansion. Hence, after Fourier transformation 
Eq.dMl) gives 



1.5 



DyO = iqXDyyD, 



where 
Dili, 



2T 



-iuj + Dq^ + Tso 



2r 



-iuj + Dq^ 



(27) 



(28) 



When expressed through Dij, the corresponding 
Cooperon components are 

Kziu^q) = -K^,{-CJ,q) - -lDyo{uJ,q) (29) 

and Kxz = K^^. At the same time 

XL(w,(Z) = KA~^,q) = -DooiLO,q) (30) 

Further, keeping only the leading terms with respect 
to small parameters fir, lot, qr and hkT, from Eqs. pop. 
(I^Tj) and one can easy calculate the factors b and c. 
The coefficients c are given by 



^yxz 



aaa arr ^ara 

yzx yzx yxz 



'^Ozz 



'^Ozz 



ttNf 
ttNf 



(31) 
(32) 



In Eq. (|31l) these coefficients change their signs with each 
permutation of lowercase indexes. The factors b, in their 
turn, are defined by 



Lrrr irrr i^aaa laaa 



' 2m*r3 



i/rrr 



-Laaa 



i{q 



a^kpTrNp 
' 4m*r4 



(33) 
(34) 
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FIG. 4: (Color online) Coordinate dependence of two phase 
shifted components of the spin-Hall current, as defined by 
Eqs. dSZl) and (|4l|. rso/ETh=0.l (red), 1 (blue),and 10 
(green). The curves are calculated at nkBT/ Em ~ 0.5 and 
eV/ETh = 1 



and the vertex functions ji calculated from Eqs. (H 
(P7)) are given by 



and 



jo(2r!,Q) 



jyi2n,Q) 



0-2 P 



1 



2m*r2 

^oo(2r!,g), 

^Dyyi2n,Q). 



2r 



Dyy{2n,Q) 



(35) 
(36) 



A following important property of Ji , J2 and 
Jiscr,>/2scr Calculated with functions and coefficients de- 
fined by Eqs. ([27| -([36 |) takes place: all these partial con- 
tributions to the spin-Hall current, after initial increasing 
with the spin-orbit coupling a, saturate when a — 00. 
This behavior is already seen in jo{2n,Q). Indeed, as 
follows from Eqs. ([55)1 and (1^51) . due to cancellation at 

— > and Q of two terms in brackets of (1551) this 



function becomes constant at large Fjo ^ a . It can 
be checked that the same combination as in brackets of 
Ea. (l35p enters into all terms contributing to the spin- 
Hall current. This sort of cancellation takes place also 
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FIG. 5: (Color online) Spin-Hall current as a function of the 
bias voltage at Tso/Eth = 1 and a; = 0. TrksT / ETh~0.25 
(red), 0.5 (blue), and 1 (green). 




x/L 

FIG. 6: (Color online) Coordinate dependence of two phase 
shifted components of the spin polarization, as defined by Eqs. 
(|39l) and rso/ETh=0.1 (red), 1 (blue), and 10 (green). 

The curves are calculated at nkBT/ Eth ~ 0.5 and eV/ Eth ~ 
1 



in normal systems and is inherent to all linear in k SOI 
couplings. In normal systems it results in the vanishing 
stationary spin-Hall effect^. Indeed, there the spin cur- 
rent is driven by the stationary electric field which, due to 
charge screening, at the small screening length is homoge- 
neous in samples of simple geometries. That guarantees 
Q = and, hence, the vanishing spin-Hall current. In 
contrast, in the considered here case of a nonstationary 
and inhomogeneous electron transport the spin-Hall cur- 
rent remains finite. The reason is that it is determined 
by the superconducting proximity effect, as it is discussed 
in the end of Section II. Consequently, a finite penetra- 
tion range of Cooper pairs results in the finite Q. More- 
over, there are also the "surface" terms, like the J2-term, 
localized near superconducting leads within the spin re- 
laxation length y^D/Tso- For these terms DQ^ ^ Tso- 
Therefore, they are not expected to saturate with larger 
SOI. 

We will normalize the spin-Hall current density accord- 
ing to 



Jc 171*0^ ^ 

T — T 



(37) 



where Jc is the critical Josephson current density defined 
by the sum over Matsubara frequencies w — TTkBT{2n+l) 



=06 



Jr = 



f 



\J± 

|A|2. 



;Doo{2ioJ,L), (38) 



The diffuson in this equation is obtained as a Fourier 
transform from and has the form Doo{2iuj, L) — 

r exp{—K,L)/Dn, with k, = \f2LjjD. The so defined 
dimensionless factor / is of the order of 1. We note 
that although a can be quite large in narrow gap 
semiconductors^^, as well as in some metallic systems^^, 
the parameter m*cP'jT in Eq. [37] is small, because the 
diffusion approximation requires hp a <C F. 
The spin density is normalized as 



Sy — SyoP , 



(39) 



where Syo is the spin polarization induced by the critical 
Josephson current in the stationary case. This polariza- 
tion is given by^ 



5. 



yo 



OLT 

2^ 



(40) 
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FIG. 7: (Color online) Spin polarization as a function of the 
bias voltage at Vso/Eth ~ 1 and a; = 0. nkBT / ETh = ^ 
(green), 0.5 (blue), and 0.1 (red). 

Since, according to (fT2|). integrand in (fT4l) contains terms 
proportional to 5n,±eV, the normalized spin current and 
spin polarization can in general be represented as 

/ = -Iism{2eVt) + l2Cos{2eVt) (41) 
P = Pi sin(2eyt) - P2 cos{2eVt) . (42) 

The spatial distribution of Ii and I2 is shown in Fig. 
4 at different values of the spin-orbit couplings. It is 
seen that the magnitude and direction of the spin-Hall 
current vary fast in the region of contacts. We recall in 
this connection that the contacts are assumed relatively 
narrow in Fig. 1 and are approximated by point-like 
sources placed at 2; = ±L/2. At some moments of time 
the normalized spin current (I42p changes its sign also as 
a function of r^o, as one can see from comparison of I2 
curves at T so = Eth and V so = WExh- Such changes are 
associated with discussed above cancellation of various 
contributions to the spin current at eV and ksT <^ T so- 
li have been pointed out that there are two types of 
terms, "bulk" and "surface" ones. At large Tgo the latter 
are localized near contacts, while the former penetrate 
deeper between and outside contacts. These qualitative 



features are clearly seen in Fig. 4. As expected, due to 
the spin-current saturation, the bulk contribution to the 
normalized current decreases with larger SOI. Indeed, Ii 
is reduced in the middle of the junction at Vso = lOpTh , 
while smaller changes are seen just near the contacts. At 
the same time, such a reduction is not so fast, as it was 
expected at first sight. At least, even at Tso = lOExh 
there is no considerable reduction of /2 at a; = 0. 

It is important to note that the spin current is not 
zero at > 1/2, while the Josephson current is absent 
there. In this region the former is driven by the time 
dependent potential, associated with the charge imbal- 
ance, rather than by a direct conversion of singlet Cooper 
pairs to triplet ones. In this spatial region the spin cur- 
rent is contributed by both "surface" and "bulk" terms. 
In Fig. 4 it extends outside the junction over the range 
^ L, because the both characteristic lengths y/D/kBT 
and y^DjeV are taken of the order of L. 

The voltage dependence of the spin-Hall current at x — 
and various temperatures is presented in Fig. 5. We 
note that both phase shifted components change sign at 
some voltages. At the same time, the behavior at the 
small bias is given by Ii ~ V'^ and /2 '-^ T^, as can also 
be checked analytically. 

The coordinate dependence of the normalized spin po- 
larization is shown in Fig. 6. In contrast to the con- 
sidered previously stationary case^, it is finite in the re- 
gion > 1/2. Similar to the spin-Hall current, such 
a behavior can be explained by the charge imbalance ef- 
fect. This effect becomes weaker at larger Tso- Indeed, at 
Tso — IO-Et/i the magnitude of P considerably increases 
inside the junction, while it becomes smaller outside it 
at l^l > 1/2. This is opposite to the spin current trend 
observed in Fig. 4. The reason is that the slow varying 
"bulk" terms are not suppressed in P at larger SOI, be- 
cause in contrast to the spin current, there is no cancel- 
lation of the spin polarization at DQ- and il <C Tso- 

A variation of the polarization with V is presented 
in Fig. 7 for Tso = Exh ■ P2 linearly turns to at 

V ^ 0, while Pi reaches its maximum there. Our nu- 
merical results show also that the magnitude of Pi at 

V = increases with Tso 00, reaching 1, that can also 
be checked analytically. At the same time, the spatial 
dependence of Pi takes the form of a step function, such 
that Pi = 1 when \x\ < 0.5L and Pi = at > 0.5L. 
Hence, as expected, in this limiting case the spin polariza- 
tion coincides with polarization calculated in the regime 
of the stationary Josephson effect^, after substitution in 
(H^ of the phase factor sm{2eVt) by sin (j), where (j) is the 
phase difference between superconductors. 

In the considered here case the spin-Hall current car- 
ries the z-oriented time dependent spin polarization in 
the transverse (y) direction. Hence, the corresponding 
time dependent spin density can be accumulated at some 
distance near fianks of the junction, while in its bulk the 
y-polarization is finite between and outside contacts. The 
spin polarization might be detected by various methods 
employed in the case of the ordinary SHE^'**. Quite ef- 
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ficient is an all-electric method used in Ref. |3- In this 
method a nonequilibrium spin polarization in the nor- 
mal metal diffuses into an adjacent ferromagnet. On the 
other hand, it is well known^**, that the spin flux through 
a ferromagnet-normal metal interface induces a voltage 
difference across the interface, that can be measured. If 
we will try, however, to extend this method to the su- 
perconducting transport, we will face a problem of eval- 
uating this voltage. It is well known how to calculate 
it in the case of a nonequilibrium flux of single particle 
spins. Much less, however, is known how to do this in 
our case, when triplet Cooper pairs contribute to this 
flux. Therefore, additional studies are necessary for a 
more complicated system than considered here. 

This work has been supported by Taiwan NSC (Con- 
tract No. 96-2112-M-009-0038-MY3) and MOE-ATU 
grant. 



Appendix A: Derivation of basic equations 

In this Abstract some explanations will be done of 
how Eqs. ([Ti| and p^ - (P5|) have been derived within 



J 



Keldysh diagrammatic technique. We start from defini- 
tions of Green functions entering into perturbation ex- 
pansions. Each of the Green functions is the matrix in 
the Keldysh spacei^: 



G = 







(Al) 



The elements of this matrix are 2x2 matrices in the 
Nambu space and 2x2 matrices in the spin space. Be- 
sides, they depend on two time arguments and two 
wavevectors k and k'. These vectors are not equal, be- 
cause there is no momentum conservation in the presence 
of impurity scattering. At the same time, the functions 
averaged over impurity positions become diagonal in the 
momentum space. Considering G as matrices in fc-space 
and also the tunneling amplitudes tk,k' in Eq. ([T]) as el- 
ements of the matrix i, one can write the expression for 
the spin current (j4]) in the fourth perturbation order with 
respect to the tunneling Hamiltonian ([T]) in the form 



Jsix,t) 



Y\ dt,Tr[a' {—Gn{t - ti)t6i2(ii, ^2)^622(^2 - t3)id2i{t3,h)iGn{t4 - t)), 



(A2) 



r 



where G and G denote the Green function of the super- 
conductor and the normal metal, respectively. The oper- 
ator k = fcj,5k,k'- We explicitly wrote the Nambu labels 
of functions, so that only the trace over spin and mo- 
mentum variables must be taken in Eq. (|A2[) . Since only 
the Josephson tunneling is considered, in the perturba- 
tion expansion we take into account only anomalous G12 

and G21 functions of superconducting leads, neglecting 
thus the usual stationary one-particle tunneling. These 
functions can be written in the form 

Gi2{t,t') ^ J-{t-t')exp{±ieV{t + t')/2) , (A3) 

where the signs "-" and relate to the left and right 

superconducting leads, respectively. G21 is obtained 
from this equation with the substitution J- J-^ and 
V ^ —V. After Fourier transform of Green functions 
in Eq. (|A2p . the time dependent exponential factors in 

Gi2(t,i') and G2i{t,t') give the factors exp(±2ieW) in 
Eq. ((131). It is important that the functions T and J^^ 
in Eq. (|A2p belong to different superconducting leads. 
Therefore, they are independently averaged over impu- 
rity positions and, hence, are diagonal in the momentum 
space. In the spin space they are proportional to the 
Pauli matrix CTz, because, as it is discussed in Section 



II, these functions are associated with the singlet Cooper 
pairing. 

It is easy to see that the Keldysh component of a prod- 
uct ABCD... of matrices having the triangular form (lAip 
can be written as the sum of products: (v4 -^5° G° £)"...) + 
{A"- CD"- ...) + {A'- B'-C'^ D" ...) + .... Ixi a\\ these prod- 
ucts the Keldysh function enters only once, while re- 
tarded and advanced functions are placed on the left and 
on the right from it, respectively. The time Fourier ex- 
pansions of thermally equilibrium Keldysh functions G 
and F can be expressed in terms of retarded and ad- 
vanced functions as 



G 



K 



(w) = (g"'(w) - G°(w)) tanh 



2fcsT 



(A4) 



We thus will apply the above expressions to the Keldysh 
component of the product in Ea. (IA2[) . This equation 
can be further simplified taking into a ccoun t that fc^T, 
eV <C A and distances of interest ^ ^JD/ IS.. Therefore, 
the main contribution to Ea. (|A2p is given by small fre- 
quencies cj ^ A. Since the difference of retarded and 
advanced anomalous functions F in Eq. (jA4p gives an 



expression proportional to 5 



[El 



A' 



j^), the corre- 



sponding Keldysh component can be neglected for small 
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The next step is averaging over disorder. From Eqs. 
(|A2|) . (|A4|) . and taking into account that G ^ az we get 
the foUowing products to be averaged G\j^azG22<^zGli, 
G1i(7zG22<^zG1i, GiiazG22<^zGii, and G\i(TzG22<7zG1i. 
In its turn, each of the Green functions is a sum of prod- 
ucts GkUk,kiGkiUk^^k2---j where Gk are diagonal in k un- 
perturbed functions and the number of impurity scatter- 
ing amphtudes U in each product is equal to the cor- 
responding perturbation order. Calculation of such av- 
erages is described in many textbooks (for example, see 
[l8| ). Briefly, within the Born approximation, assuming 
random impurity positions the averages of Uk,k, products 
decouples into pair averages. In this way each pair enters 
as an effective two-particle interaction carrying a zero fre- 
quency. So, the average of a Green function can be ex- 
pressed through the self-energy. In Eq. |6]this self-energy 
is given by iV. When a product of Green functions is 
averaged, a considerable simplification takes place in the 
semiclassical approximation when T <^ Ep. In this case 
a special class of the so called "ladder" diagrams domi- 
nates in the perturbation expansion over disorder. Let us 
consider, for example the average {G'[iazG22crzGii)iaip- 
Combining in pairs U-s in G22 with its neighbors on the 



left and on the right we obtain two ladder series. There 
are "Cooperons", according to definition ([9|). They are 
shown as gray boxes in Fig. 2a. At the same time, one 
can not build the ladder out of the pair G^i Gi^. The 
reason is that the sum over k of a typical ladder element, 
the product G^\^ (g) G^^^, where G°'' is given by Eq. ©, 
turns to zero, because both functions have poles in the 
same semiplane of the complex variable Ek ■ On the same 
reason the ladders built of G22 ^ Gfi also turn to zero, 
as follows from definition ([6|). Therefore, the average 
{GliazG 22!^ zGfi) imp contains only one Cooperon origi- 
nating from the first two functions. Besides, the combi- 
nation G\^ Gil entering into this product, results in a 
diffuson defined by Eq. ([8]). The corresponding Feynman 
diagram is shown in Fig. 2b. 

In the same way, as the spin current, one may calculate 
the electric charge, substituting in Eq. (jA2|) a^ky/2m* 
for e. This charge is given by the sum of polygons in 
Figs. 3a and 3b. They represent Jiscr and J2scr, respec- 
tively. Further, the screening electric potential is cal- 
culated within the random phase approximation. This 
potential drives the spin-Hall effect in the same way as 
in normal systems.^ 
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